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I discuss isospin violation in an effective field theory appraoch of the nucleon- 
nucleon interaction. The observed leading order charge independence breaking 
is explained in terms of one— pion exchange together with a four-nucleon contact 
term. A general classification scheme of corrections to this leading order results 
is presented and an explanation is offered why charge-dependent pion— nucleon 
coupling constants can not play a role in the observed isospin violation in the NN 
scattering lengths. Charge symmetry breaking, the size of ^-exchange graphs 
and the pionless theory are also discussed. 



1 Introduction 

Heisenberg introduced the notion of isospin in 1930 as a label to distinguish 
between neutrons and protons when he was studying atomic nuclei. He noticed 
that the clustering of stable nuclei around N ~ Z implied a new (internal) 
symmetry. This was later extended to the pion triplet, which is believed to 
mediate the long range part of the nuclear force. Heisenberg's ideas were 
used by Breit and collaborators (and others) to establish the so-called charge 
independence, observing that after removal of the Coulomb effects between 
protons, the pp, np and nn force were equal to a good precision. A weaker 
condition is the so-called charge symmetry, which states that the pp and nn 
force should be equal with no assertion to the np force. Since that time, 
QCD has emerged as the theory of the strong interactions. In the absence of 
electroweak interactions, isospin symmetry is observed in the sector of the two 
lightest quark flavors for equal masses, m u = m^, i.e. the QCD Hamiltonian is 
invariant under a global (flavor) transformation of the type 

fj^q' = Vq = V (f) , VeSU(2) . (1) 

Of course, in nature the light quark masses are not equal. In fact, they differ 
considerably, n% u /m,d = 0.55. Still, isospin is a good approximate symmetry 
because what counts is not the ratio of the quark masses but rather their differ- 
ence compared to the typical scale of strong interactions, (m^ — to u )/Aqcd 
1. QCD thus provides a reason why isospin is such a good approximate sym- 
metry. This can hardly be considered "accidental" . In the presence of electro- 
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magnetism, further isospin violation is induced since the charges of the quarks 
are unequal. These effects arc in general small due to the explicit apperance of 
the fine structure constant a ~ 1/137. In fact, in many cases the strong and 
electromagnetic isospin-violating effects are of the same size. A good example 
is the neutron-proton mass difference. On the quark level, charge symmetry 
is a rotation about the 2-axis in isospin space 

"-OH" 2 '- C0=(~«) • pi=i ■ (2) 

Our aim is get some insight about isospin symmetry (charge independence) 
and charge symmetry in the two-nucleon system based on an effective field 
theory approach. On the level of the two-nucleon force, charge independence 
implies a potential of the form 

Vij = A + Bt ({) ■ T(j) , (3) 

where A and B are functions of the nucleon spin and momentum (or space) 
coordinates and '(i),(J)' labels the interacting nucleons. Charge symmetry 
allows for a more general two-body force, 

% cs = A + Bt (i) ■ t U} + Crf t) ■ rf f) . (4) 

If one works out the Coulomb potential in terms of the nucleon charge matrix 
Q = e(l + t 3 )/2, it is obvious the electromagnetic effects lead to breaking of 
charge independence and charge symmetry, 

e 2 

^.Coulomb = ^ 0- + T h ) (! + T U) ) ■ ( 5 ) 

I now wish to show how the question of isospin violation can be addressed 
in an effective field theory approach to the NN interactions. There are twa 
equivalent ways to proceed. One is to follow the suggestion of Weinberg El 
and perform the chiral expansion on the level of the potential (i.e. for the 
irreducible diagrams) and iterate this potential in the Lippmann-Schwinger 
equation to generate the bound and scattering states. The other is to directly 
expand the amplitude after performing a nonperturbativejresummation for the 
S-wavesEl The first path has been followed by van KolckB and I will consider 
the second option hereQ Pions are included and treated perturbatively in this 
approach. For the topic addressed here, the discussion whether or not pions 
should be treated in such a way is inessential. In fact, as shown in section 6, 
one can even work with a pionless theory. In that case, however, the physical 
interpretation of CIB and CSB becomes much less transparent. 
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2 Facts, explanations and mysteries 

It is well established that_the nucleon-nucleon interactions are charge depen- 
dent (for a review, see e.g0 and a recent summary is given inLT). For example, 
in the 1 Sq channel one has for the scattering lengths a and the effective ranges 
r 



Aa C iB = - (a™ + a pp ) - a np = 5.7 ± 0.3 fm , (6) 
Ar C iB = i (r nn + r pp ) - r np = 0.05 ± 0.08 fm . (7) 



These numbers for charge independence breaking (CIB) are based on the Ni- 
jmegen potential!] and the Coulomb effects for pp (nn) scattering are subtracted 
basecLon standard methods (for a treatment of this in an EFT framework see 
refslfm). The charge independence breaking in the scattering lengths is large, 
of the order of 25%, since a np = (—23.714 ± 0.013) fm. Of course, it is mag- 
nified at threshold due to kinematic factors (as witnessed by the disapperance 
of the effect in the effective range). In addition, there are charge symmetry 
breaking (CSB) effects leading to different values for the pp and nn threshold 
parameters, 

AacsB = a pp -a nn = 1.5±0.5 fm , ArcsB = r pp — r nn = 0.10±0.12 fm . (8) 

Combining these numbers gives as central values a nn = — 18.8fm and a pp = 
— 17.3fm. Both the CIB and CSB effects have been studied intensively within 
potential models of the nucleon-nucleon (NN) interactions. In such approaches, 
the dominant CIB comes from the charged to neutral pion mass difference in 
the one-pion exchange (OPE), Aag^ E ~ 3.6 ± 0.2 fm. Additional contribu- 
tions come from jit and 2ir (TPE) exchanges. According to some calculations, 
their size is approximately l/3 rd of the leading OPE contribution. For the 
TPE, significantly smaller results can also be found in the literature— see e.g. 
table 3.3 in rem In case of the nj contribution, a recent calculation!^ gives a 
value which is more than a factor of ten smaller than the leading OPE contri- 
bution. We will come back to these issues later (see section 5). Note also that 
the charge dependence in the pion-nucleon coupling constants (if existing) in 
OPE and TPE almost entirely cancel. Naively, one would expect a possible 
charge dependence of the pion-nucleon coupling constants to be of a similar 
importance as the pion mass difference in OPE. This strong suppression of 
charge-dependent couplings has so far been eluded a deeper understanding. 
We will come back to this later on. In the meson-exchange picture, CSB orig- 
inates mostly from p — lu mixing, Aa^gg ~ 1.2 ± 0.4 fm. Other contributions 
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due to 7T — rj, n—r]' mixing or the proton-neutron mass difference are known to 
be much smaller. It will be shown in what follows that EFT gives novel insight 
concerning the size of the TPE and 717 contributions as well as the suppression 
of possible charge-dependent nN coupling constants. 

3 Effective Lagrangian and power counting 

The effective Lagrangian underlying the analysis of isospin violation in the 
two-nucleon system consists of a string of terms, 

Lett — Ctttt + C-nN + £-NN ■ (9) 

Each of these terms splits into a strong and an electromagnetic one. The 
strong contributions include (besides many others) the isospin-violation due 
to the quark mass difference m M — whereas the em terms stem from the 
dynamics of the virtual photons (photon loops and em counterterms). To 
include virtual photons-Ln the pion and the pion-nucleon system is by now a 
standard procedure Onia We note that in the pion and pion-nucleon sector, 
one can effectively count the electric charge as a small momentum or meson 
mass. This is based on the observation that M„/A x ~ e/\/47r = y/a ~ 1/10 
since A x ~ 4-71"/^ = 1.2 GeV. It is thus possible to turn the dual expansion in 
small momenta/meson masses on one side and in the electric coupling e on 
the other side into an expansion with one generic small parameter. We also 
remark that from here we use the fine structure constant a = e 2 /Air as the em 
expansion parameter. Because of the additional 1/47T factor (from integrating 
out hard photons) one expects the corresponding em constants to be of order 
one (of natural size) dimensional analysis of such terms in the.pion-nucleon 
sector is given in refEj and similarly for the pion sector in refsLUIiil. 

We now turn to the two-nucleon sector. One can either expand the poten- 
tial (irreducible diagrams) or directly the amplitudes. A complication arises 
due to the unnaturally large S-wave scattering lengths, 1/a <C M n . A sys- 
tematic power counting scheme to deal with this directly on the level of the 
scattering amplitude is the recently proposed power-divergence subtraction 
method (PDS) of Kaplan, Savage and Wise (KSW) H.g Essentially, one re- 
sums the lowest order local four-nucleon contact terms ~ Cq(N'N) 2 (in the 
S-waves) to generate the large scattering lengths and treats the remaining ef- 
fects perturbatively, in particular also pion exchange. This means that most 

a There exist by now modifications of this approach and it as been argued that it is equivalent 
to cut— off schemes. We do not want to enter this discussion here but rather stick its original 
version. 
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low-energy observables are dominated by contact interactions. The chiral ex- 
pansion for NN scattering entails a new scale n of the order of 300 MeV, 
so that one can systematically treat external momenta up to the size of the 
pion mass. There have_been suggestions that the radius of convergence can 
be somewhat enlarged E3, but in any case Ann is considerably smaller than 
the typical scale of about 1 GeV appearing in the pion-nucleon sector. It is 
straigthforward to extend the underlying effective Lagrangian to include the 
strong isospin- violating and electromagnetic four-fermion contact interactions. 
Consider first the strong terms. Up to one quark mass insertion, the effective 
Lagrangian takes the form (note that terms involving the Pauli isospin matrices 
t can be obtained by Fierz reordering) 

£f N = h(N^N) 2 + l 2 {N^aNf + l 3 (N^ (x+)N) (N^N) + h(N^ x+N)(N* N) 
+ l 5 (N^a( X+ )N)(N^aN) + l 6 (N^a X +N)(N^aN) + ... , (10) 

where (. . .) denotes the trace in flavor space and \+ contains the light quark 
masses. \+ = X+ — (x+)/2 is only non-vanishing for m u ^ rrid and the ellipsis 
denotes terms with two (or more) derivatives acting on the nucleon fields. 
Similarly, one can construct the em terms. The ones without derivatives on 
the nucleon fields read (possible terms involving the Pauli isospin matrices are 
of no relevance for the following) 

£nn = & {n(Q 2 + - Q 2 _) + r 2 Q + (Q+)} N(N^N) 
+ rfcr {r 3 (Q 2 + - Q 2 _) + r 4 Q+(Q+)} N(N^aN) 

+ AT f {r 5 Q+ + r 6 (Q+)}N(N^Q + N) 

+ N^{r 7 Q + + r s {Q+)}N(N^aQ + N) 

+ r g (N^Q + N) 2 + r w {N^aQ + N) 2 , (11) 

with Q± = uQ^u ± u'Qu' and Q the nucleon charge matrix. All these terms 
are due to the integration of the "hard" virtual photons. The "soft" ones are 
contained in the periipent covariant derivatives acting on the pion or nucleon 
fields, see e.g. refsll3li3. We remark that since Ann is significantly smaller than 
A x , it does not pay to treat the expansion in the generic KSW momentum^] 

Q/A NN ~ ft/Ajvjv ~ M n /A NN (12) 

'So that no confusion with the charge matrix can arise, we denote this parameter with an 
extra "tilde" contrary to common use. A pedagogical discussion of this expansion parameter 
can be found in refta. 
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(with u the scale of dimensional regularization) simultaneously with the one 
in the fine structure constant (as it is done e.g. in the pion-nucleon sector). 
Instead, one has to assign to each term a double expansion parameter, 

{Q/AnnY 1 a rn , Q/ Ann =s 1/3 , a = 1/137 , (13) 

with n and m integers (m > 1 and n can be negative). Clearly, even if we 
work at next-to-next-to leading order (NNLO) which entails a supression by 
two powers in Q, we need only be concerned with the leading em terms ~ a 
since a <C 1. Lowest order charge independence breaking is due to a term 
~ (N^t 3 N) 2 whereas charge symmetry breaking at that order is given by a 
structure ~ (N't 3 N)(N 1f N). In the KSW approach, it is customary to project 
the Lagrangian terms on the pertinent NN partial waves. Denoting by (3 the 
So partial wave for a given cms energy E cma , the Born amplitudes for the 
lowest order CIB and CSB operators between the various two-nucleon states 
takes the form 



(/3,pp|£^|/3,pp) = - 
(/3,nn\C e ^ N \(3,nn) =- 




where the coupling constant Eq ' (/i) will be determined later on. From the 
renormalization group equations (RGE) it will also be possible to derive the 

( 1 2) 

scaling properties of the E 2r ^ ' (//) for four-nucleon contact operators with 2n 
derivatives. The terms with the superscript '(1)' refer to CIB whereas the 
second ones relevant for CSB are denoted by the superscript '(2)'- Higher 
order operators are denoted accordingly. There is, of course, also a CIB con- 
tribution to the np matrix element. To be consistent with the charge sym- 
metric calculation of rem, one has to absorb its effect in the constant D2, 
i.e. it amounts to a finite renormalization of D2 and is thus not observable. 
In eq.(|lj), p — V ME cms is the nucleon cms momentum and M the nucleon 
mass. 

4 Isospin violation in nucleon nucleon scattering 

Let me now work out-the pertinent 1 5'o phase shifts. The np case has already 
been calculated in refH, the leading amplitude being of order Q~ l . We want to 
calculate the leading em corrections (here, "Cs" means Coulomb-subtracted 
and we neglect the small magnetic corrections for the nn system), 
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Figure 1: Relevant graphs contributing to charge independence breaking at leading order 
aQ~ 2 . The blob stands for the resummation of the lowest order (N> TV) 2 contact term ~ Co- 
The open (filled) circle denotes a pion mass insertion ~ Sm 2 (an insertion of the leading 

four-nucleon operators ~ aE^). For charge symmetry breaking, the leading contribution 

(2) 

is given by the last diagram with the filled circle denoting an insertion ~ ctEg . 

The pertinent diagrams are shown in fig. 1. Besides the contact terms written 
down in eq. (|l4|) , we have to consider the pion mass difference in one-pion 
exchange. OPE diagrams with different pion masses have the isospin structure 
012 = T (i)a^x br (i)6 = (a + &)t(i) f( 2 ) — r 3 2 . and obviously lead to CIB. 
From the first graph in fig. 1 one easily concludes that these effects are of order 
aQ~ 2 since the diagram contains two derivative irN couplings of order Q and 
two pion propagators of order Q~ 2 . The notation follows closely the one of 
KSW in that the corresponding three amplitudes are called At^^ 1 '™ ■ Here, 
the first (second) subscript refers to the power in a (Q) and the superscripts 
to the first three diagrams of the figure. This gives 
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with A-i 



Ao,-i the leading term in the expansion of the np 1 So amplitude 

Co 



(17) 



Furthermore, /„■ — 92.4 MeV is the pion decay constant, qa — 1-26 the axial- 
vector coupling constant and m 2 = 2m 2 ± — rn^ the "average" pion mass. 
The somewhat unusually looking expressions for m 2 and 8m 2 are derived 
in appendix A. Note that while the diagrams II and III are convergent, IV 
diverges logarithmically. Therefore, the Lagrangian must contain a coun- 
terterm of the structure (p)a{N Ji t 3 N){N^ t s N) to make the amplitude 
scale-independent. The leading CIB contribution thus scales as Q~ 2 and for 
higher order operators with In derivatives we can establish the scaling prop- 
erty i?2n ~ Q~ 2+n . This does not contradict the KSW power counting for the 
isospin symmetric theory since aCl. The insertion from this contact term is 
shown in the last diagram of fig. 1 and leads to an additional contribution to 
AA. In complete analogy, we can treat the leading order CSB effect which is 
due to an operator of the form aE^ (N Jt T3N)(N^ N). This term is, however, 
finite. Putting pieces together, this gives for the nn and pp case 
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where the coupling constants £q ' (p.) obey the renormalization group equa- 
tions, 



pa 
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dE { j\p) M (1) 



dp 
dE<£\p) 



^E^Ui)C (p)p-Sm 2 ^^C 2 (p) . 
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(20) 



The relation of the pp and nn scattering lengths to the np one reads (of course, 
in the pp system Coulomb subtraction is assumed), 



1_ = J_ _ AnajE^' +E^>) 
- ~ a n/n MCI 

(i) _ P (2)> 



1 1 47ra(^ 1J - S^) 



MC 2 



+ A 



(21) 
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2 g\(-C M(m - 2u) + 8tt + CpMrn ln(r7t 2 / /i 2 )) 



A = 5m ^ : jr-Pi 

12mnC f4 



For the effective ranges, we have oniy CIB 



,g 2 A (C Mn + 47r)(C M(3 J u - 2m) + 12tt) 
67tMto 4 C2/2 



Note that the relations in eqs.( |2l|j22j ) are perturbatively scale-independent and 
that the last one does not contain any unknown parameter. We remark that for 
the CIB scattering lengths difference the pion contribution alone is not scale- 
independent and can thus never be uniquely disentangled from the contact 
term contribution ~ Eq . While the leading OPE contribution resembles the 
result obtained in meson exchange models, the mandatory appearance of this 
contact term is a distinctively new feature of the effective field theory approach. 
Matters are somewhat different for CSB. Here, to leading order there is simply 

(2) 

a four-nucleon contact term proportional to the constant E . Its value can 
be determined from a fit to the empirical value given in cq.(|^). Corrections 
to the leading order CIB and CSB effect are classified below. We add the 
following remarks. Eqs.(^l|) and (p2|) are not exact but hold perturbatively. 
Similarly, the phase shifts are also defined perturbatively and therefore one 
can not have an exact effective range expansion of 8 = So + #i (referring to the 
isospin symmetric case for the moment). If one works e.g. at next -to-leading 
order, the norm of the S-matrix is not exactly one, S = 1+0(Q 2 ) ^ 1. From 
this it follows that if one adjusts the parameters so as to reproduce the phase 
shifts in a certain momentum range, one can only approximately reproduce 
the scattering length and the effective range, a ~ a oxp and r ~ r cxp , although 
the phase shifts may be very well described. This should be kept in mind. 

We are now in the position to calculate the phase shifts. In total, we have 
five coupling constants, three of them related to the isospin-symmetric case, 
called Co, C 2 and D 2 . We fit the paramters Gq^Di to the np phase shift 
under the condition that the scattering length and effective range are exactly 
reproduced. The two new parameters Eq ' are determined from the nn and 
pp scattering lengths. The resulting parameters at p — m are of natural size, 

C = -3.46 fm 2 , C 2 = 2.75 fm 4 , D 2 = 0.07 fm 4 , 
Ef> = -6.47 fm 2 , E (2) = 1.10 fm 2 . (23) 

The np and nn 1 Sq phase shifts are shown in fig. 2. To arrive at these curve, the 

physical masses for the proton and the neutron have been used. We stress again 

fi 2) 

that the effect of the em terms ~ Eq is small because of the explicit factor of 
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a not shown in eg. (p3|) . Having fixed these parameters, we can now predict the 
nn 1 Sq phase shift as depicted by the solid line in fig. 2. It agress with nn phase 
shift extracted from the Argonne V18 potential (with the scattering length and 
effective range exactly reproduced) up to momemta of about 100 MeV. The 
analogous curve for the pp system (not shown in the figure) is close to the solid 
line since the CSB effects are very small. 
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Figure 2: 1 So phase shifts for the np (dashed line) and nn (solid line) systems versus the 
nucleon cms momentum. The emipirical values for the np case (open squares) are taken 
from the Nijmegen analysis. The open octagons are the nn "data" based on the Argonne 
V18 potential. 



5 Classification scheme 

In the framework presented here, it is straightforward to work out the vari- 
ous leading (LO), next-to-leading (NLO) and next-to-next-to-leading order 
(NNLO) contributions to CIB and CSB, with respect to the expansion in Q, 
to leading order in a and the light quark mass difference. Here, I will simply 
enumerate the pertinent contributions which appear at a given order. This nat- 
urally also gives an estimate about their relative numerical importance based 
on simple dimensional analysis. Based on such an analysis, one can shed some 
light on the puzzles discussed in section 2, in particular the size of the ttj 
graphs and the role of charge-dependent nN coupling constants. Let me first 
consider CIB. In the classification scheme given below, TPE/3PE stands for 
two /three-pion-exchange. 
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LO 

aQ- 2 


Pion mass difference in OPE, 

four-nucleon contact interaction with no derivatives. 


NLO 

aQ~ x , eQ- 1 


Pion mass difference in TPE, 

four-nucleon contact interaction with two derivatives, 
insertions proportional to the np mass difference, 


NNLO 


Pion mass difference in 3PE, 

four-nucleon contact interaction with four derivatives, 
CIB in the pion-nucleon coupling constants, 
7T7-exchange. 



Note that there are no CIB effects due to the light quark mass difference linear 
in e = m u — rrid from OPE and from the four-nucleon contact terms. These 
type of terms can only appear to second order in e. Effects due to the charge 
dependence of the pion-nucleon coupling constants, i.e. isospin breaking terms 
from C^, only start to contribute at order aQ . Such effects are therefore 
suppressed by two orders of Q compared to the leading terms, i.e. by a nu- 
merical factor of about (1/3) 2 . This finding is in agreement with the various 
numerical analyses performed in potential models (if one assumes that there 
is indeed some charge dependence in the pion-nucleon coupling constants, see 
e.g. the discussion in reffl.) We remark that the TPE contribution is nomi- 
nally suppressed by one order in the small parameter, i.e. one would expect 
its contribution to be one third of the leading OPE with different pion masses. 
This agrees with some (but not all) meson-exchange model calculations. 



a b c 



Figure 3: 7T7 graphs. For CIB, the leading contribution comes from the three diagrams a), 
b) and c). Due to the Kroll-Ruderman vertex in a) and c) the exchanged pion is charged. 
In case of CSB, the dominant contribution comes from b) with neutral pion exchange, c) 
vanishes due to parity and a) is suppressed by two additional powers in Q. Crossed graphs 
are not shown. Also omitted are the initial and final-state interactions between the nucleon. 

Some remarks on the 77r graphs (some of them shown in fig. 3) are in order. 
First, we note that their contribution for CIB is suppressed by two powers 
in Q with respect to the leading terms (the power counting for these graphs 
is detailed in appendix B). The larger contribution found in some meson- 
exchange model calculations (approximately l/3 rd of the leading OPE con- 
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tribution) is therefore at first sight at variance with n dimensional analysis 
since one would expect a suppression of about (1/3) 2 U Howejier, the much 
smaller value obtained in the more recent calculation of reflij (performed 
in the static limit for the nucleons and Coulomb gauge for the photons), 
AaciB(i"7)/AaciB(OPE) ~ —1/14.6, is in much better agreement with the 
suppression by two powers in Q. This deserves further study. Second, the jtt 
graphs evaluated in heavy baryon chiraLperturbation theory have been used 
in a refined potential model calculation^ The np low-energy parameters 
were found to be very little affected. It is also instructive to see how the cor- 
rections due to the np mass difference come about. For that, let us write the 
nucleon mass as M + 5M, where SM subsumes the em and strong contribution 
to the np mass splitting, i.e. these terms are 0(a) and 0(e). Differentiating 
the leading isospin-symmertic amplitude eq.(|l7|) with respect to the nucleon 
mass shift gives 

W = -^ + ip) ' (24) 

from which it follows that 

A 2 1 - 
SA-i = -5M-=±(u + ip) ~eQ- 1 ,oQ -1 . (25) 

We now turn to CSB. The pattern of the various contributions looks different: 



LO 

aQ- 2 ,eQ- 2 


Electromagnetic four-nucleon contact interactions 
with no derivatives. 


NLO 


Em four-nucleon contact terms with two derivatives, 
strong isospin-breaking contact terms w/ two derivatives, 
insertions proportional to the np mass difference. 



We remark that the leading order CSB effects do not modify the effective range 
but the corresponding scattering length. Also, the leading ir^ graphs (with the 
exchanged pion now being neutral) are of order aQ°, i.e. suppressed by two 
orders in the KSW expansion parameter. This can be used as an argument to 
simply ignore such contributions in the study of CSB. Within meson-exchange 
models, their contribution comes out small but there seems to be no consensus 
about the actual sizeQ A similar hierarchy can also be derived based on t 
Weinberg power counting (which applies to the potential) as discussed in re: 

6 The pionless theory 

It is also of interest to consider a simplified theory in which even the pions are 
integrated out and one only deals with four-nucleon contact terms. Such an 
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approximation is clearly only justifiable in a KSW-like approach as followed 
here. We do not want to enter the discussion about how to treat the pions 
(perturbative versus nonperturbative) but rather wish to see how accurate the 
pionless theory works for the question of isospin violation. In such an approach 
the connection to the chiral aspects of QCD is, of course, lost. Integrating out 
the pions modifies the coupling constants of the various contact terms like 

EFT with pions: Co, C 2 , D 2 and , 
EFT without pions: C Q , C 2 and . 

In the pionless theory, there is no more D 2 term. One has thus four parameters 
which can be fixed from the three scattering lengths and the np effective range. 
The values of the parameters are 

Co = -3.42 fm 2 , C 2 = 6.05 fm 4 , 
E ( 1] = 8.09 fm 2 , E^ = 1.07 fm 2 . (26) 




p [MeV] 

Figure 4: x 5o phase shifts for the np (dashed line) and nn (solid line) systems versus the 
nucleon cms momentum in the pionless theory. The emipirical values for the np case (open 
squares) are taken from the Nijmegen analysis. 

Compared to the values for the theory with pions given in eq. (p3|) , we note 
that the value of Co almost the same as the one of Co . This is expected since 
we are close to a fix point. Noticeably, there is a dramatic difference in Eq 1 ^ 
compared to Eq X \ which is also expected since the former pion contribution 
is now completely absorbed in the contact term. The resulting np and nn 
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1 S'o phase shifts shown in fig. 4. The description of the phase shifts has clearly 
worsened, they start to deviate from the data for momenta larger than 50 MeV 
and are considerably off the data for higher momenta, say about 200 MeV. 
These results are not surprising, but it is comforting to see that the pionless 
theory indeed describes the small momentum region correctly. For example, 
to this order the nn effective range equals the one of the np system which is in 
agreement with the data (partial wave analysis). However, the interpretation of 
the CIB is now very different since it is entirely given by the isospin-violating 
contact interaction. This makes a comparison with the results obtained in 
potential models very difficult, if not impossible. 

7 Summary 

We have shown that isospin violation can be systematically included in the 
effective field theory approach to the two-nucleon system in the KSW formu- 
lation. For that, one has to construct the most general effective Lagrangian 
containing virtual photons and extend the power counting accordingly. We 
have shown that this framework allows one to systematically classify the vari- 
ous contributions to charge indepedence and charge symmetry breaking (CIB 
and CSB). In particular, the power counting combined with dimensional anal- 
ysis allows one to understand the suppression of contributions from a pos- 
sible charge-dependence in the pion-nucleon coupling constants. Including 
the pions, the leading CIB breaking effects are the pion mass difference in 
OPE together with a four-nucleon contact term. These effects scale as aQ~ 2 . 
Power counting lets one expect that the much debated contributions from 
two-pion exchange and graphs are suppressed by factors of 1/3 and (1/3) 2 , 
respectively. This is in agreement with some, but not all, previous more model- 
dependent calculations. The leading charge symmetry breaking is simply given 
by a four-nucleon contact term. One can also use the pionless theory to de- 
scribe the observed CIB and CSB, but the physical interpretation is much 
less clear. It would certainly be interesting to investigate such effects in the 
presence of external sources. 
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Appendix A: The average pion mass 



Here, I briefly show how the average pion mass used in the np system in 
section 4 comes about. Consider the graphs shown in fig. 5. The strong 
interaction potential for OPE in the nn and pp system is given by the standard 
Yukawa form 



Vpp — 



n 2 

9ttN 



• V)(a (2) -V)e- m « r 



47r 4M 2 r 

where m n is a shorthand notation for the neutral pion mass. 



(27) 



n 



n 













71° 











Figure 5: OPE graphs for the nn (left) and the np system (middle and right). 

For the np case, we also have charged pion exchange, 

V np = V nn [- 1 + 2e~ 5m r ]~y M [l- 25m r] , (28) 

with 5m = m± — m and m± denotes mass of the charged pion. So effectively 
one has for the np OPE a mass 

m = m + 28m or m 2 — m\ + 25m 2 , (29) 

since m 2 ~ TOq + ASmmo = TOq + 4(m± — mo)mo ~ m§ + 25m 2 (note that 
5m 2 ^ (5m) 2 ). These are exactly the formulae used in section 4. 

Appendix B: Power counting for the ttj graphs 

In this appendix, we will work out the explicit power counting of the various 
7T7-exchange diagrams (cf. fig. 3). First we need the counting rules for the 
propagators, vertices etc.: Pion propagator ~ Q~ 2 , nucleon propagator ~ 
Q^ 1 , pion-nucleon vertex ~ Q, photon-nucleon vertex ~ Q°, charged pion- 
photon-nucleon vertex ~ Q°, neutral pion-photon-nucleon vertex ~ Q and 
loop integration ~ Q 4 . This gives (displayed is the respective power in Q for 
all vertices, propagators, etc for the specified diagram): 
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Graph Vertices Propagators Loops Total 



3a (tt*) 
3a (tt ) 
3b (7r±) 
3b (tt ) 
3c (tt ± ) 
3c (tt ) 



For CIB, we have charged and neutral pion exchange and all diagrams 3a,b,c 
start contributing Q°. In the case of CSB, matters are somewhat different. 
Graph 3b is leading at Q°, while 3c is formally of order Q but vanishes due to 
parity. The contribution of graph 3a to CSB is suppressed by two additional 
powers of the small momentum. 
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